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Abstract. We present an asteroseismic study of the fast rotating star HD 187642 (Altair), recently discovered to be a d Scuti 
pulsator. We have computed models taking into account rotation for increasing rotational velocities. We investigate the relation 
between the fundamental radial mode and the first overtone in the framework of Petersen diagrams. The effects of rotation 
on such diagrams, which become important at rotational velocities above 150kms“', as well as the domain of validity of our 
seismic tools are discussed. We also investigate the radial and non-radial modes in order to constrain models fitting the five 
most dominant observed oscillation modes. 
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1. Introduction 

The bright star A ltair (a Aql) was recently observed by 
iBuzasi et al 1 j2005l) (hereafter paper-I) with the star tracker on 
the Wide-Field Infrared Explorer (WIRE) satellite. The overall 
observations span from 18 October until 12 November 1999, 
with exposures of 0.5 seconds taken during around 40% of the 
spacecraft orbital period (96 minutes). The analysis of the ob¬ 
servations made in paper-I reveals Altair to be a low-amplitude 
variable star (Am < 1 ppt), pulsating with at least 7 oscilla¬ 
tion modes. These results suggest that many other non-variable 
stars may indeed turn out to be variable when investigated with 
accurate space observations. 

Since Altair lies toward the low-mass end of the instability 
strip and no abundance anomalies or Pop II characteristics are 
shown, the authors identified it as a d Scuti star. The b Scuti 
stars are representative of intermediate mass stars with spec¬ 
tral types from A to E. They are located on and just off the 
main sequence, in the faint part of the Cepheid instability strip 
(luminosity classes V & IV). Hydrodynamicalprocesses occur¬ 
ring in stellar interiors remain poorly understood, b Scuti stars 
seem particularly suitable for the study of such physical pro¬ 
cess, eg. (a) convective overshoot from the core, which causes 
extension of the mixed region beyond the edge of the core as 
defined by the Schwarzschild criterion, affecting evolution; and 


Send offprint requests to: J.C. SuMez 
* also associated researcher at LESIA, Observatoire de Paris- 
Meudon, UMR8109 


(b) the balance between meridional circulation and rotationally 
induced turbulence gen erates chem ical mixing and angular mo¬ 
mentum redistribution (IZahrJl 992 h . 


From the observational side, great efforts have been made 
within last decades in developing the seismolo gy of S Scuti 
stars within co ordinated networks, e.g.; STE PHI jMichel et alJ 
I 2 OOOI) or DSN (lBregeil200(i lHandleii2000l) . However, several 
aspects of the pulsating behaviour of these sta rs are not com¬ 
pletely understood (see Templeton et For more de¬ 

tails, an interesting revie w of unsolv ed problems in stellar pul¬ 
sation physics is given in ICoxI 1I2OO2I1 . Due to the complexity of 
the oscillation spectra of 6 Scuti stars, the identification of de¬ 
tected modes is often difficult and require additional informa - 

reger et al.lll( 


tion (see for instance fviskiim et al.l 
A unique mode identification is often impossible and this ham¬ 
pers the seismology studies for these stars. Additional uncer¬ 
tainties arise from the effect of rapid rotation, both directly, 
on the hydrostatic balance in the star and, perhaps more im¬ 
portantly, through mixing caused by circulation or instabilities 
induced by rotation. 


Intermediate mass stars, like A type stars (5 Scuti stars, 
y Dor, etc) are known to be rapid rotators. Stars with 100 < 
vsinf < 200 km s ' are no longer spherically symmetric but are 
oblate spheroids due to the centrifugal force. Rotation modi¬ 
fies the structure of the star and thereby the propagation cavity 
of the modes. The characteristic pattern of symmetric multi- 
plets split by rotation is thus broken. In the framework of a 
linear perturbation analysis, the second order effects induce 
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Fig. 1. HR diagram containing the whole sample of models considered in this work (filled circles) as well as the observational 
photometric error box. Dotted and dot-dashed lines correspond to non-rotating evolutionary tracks of 1.70 and 1.76 Mq models 
respectively. Continuous curves represent those of 1.70Mq models with rotational velocities, from top to bottom, of 50, 100 and 
150kms“' respectively. Finally, dashed lines represent evolutionary tracks of 1.75 Mq models with rotational velocities varying, 
from top to bottom, of 50, 100 and 200 km s ' respectively. 


Table 1. Altair observed frequencies (from paper I). From left 
to right, columns provide frequencies in cd ', in /rHz, ampli¬ 
tudes in ppm, and finally, /i //j= 2,...,7 frequency ratios. 



(c/d) 

(pHz) 

A (ppm) 


/l 

15.768 

182.50 

413 


/2 

20.785 

240.56 

373 

0.759 

/a 

25.952 

300.37 

244 

0.607 

/4 

15.990 

185.07 

225 

0.986 

fs 

16.182 

187.29 

139 

0.974 

h 

23.279 

269.43 

111 

0.677 

fl 

28.408 

328.80 

132 

0.555 


strong asymmetries in the sp litting of multiplets dSai 3 ll98lt 


Dziembowski & Goodelll992 

) and shifts which cannot be ne- 

glected even for radial modes 

Soufi et al.|(|l995|). The star stud- 


ied here is a very rapidly rotating A type star. Therefore ro¬ 
tation must be taken into account, not only when computing 
equilibrium models but also in the computation of the oscil- 
lation frequencies. The forthcoming space mission COROT 
( iBaglin et alJ I 2 OO 2 I) . represents a very good opportunity for 
investigating such stars since such high frequency resolution 
data will allow us to test theoretically predicted effects of rota¬ 


tion. The paper is structured as follows: In Sect. 12 fundamen¬ 
tal parameters necessary for the modelling of Altair are given. 
Equilibrium models as well as our adiabatic oscillation code 
are described in Sect.|2 Section|4]presents a discussion of the 
two seismic approaches followed: 1) considering radial modes 
and analysing the Petersen diagrams, and 2) considering radial 
and non radial modes. We also discuss a possible modal iden¬ 
tification. Finally, major problems encountered for modelling 
Altair and conclusions are presented in Sect.|3 

2. Fundamental parameters 

Several values of the effective temperature and surface grav- 
itv of Altair (a Aql) can be found in the literature. Recently, 
tersnamer & Northll2003l) proposed Teg = 7550/T and logg = 
4.13 + 0.3 dex, derived from photometric measurements 
(Geneva system). From Hipparcos measurements (a parallax 
of 194.44 + 0.94 mas) combined with the observed V ma gni- 
tude and the bolometric correction given bv iFlowe 3 (ll996l) . we 
obtain a bolometric magnitude of Mboi =2.18+0.1 dex. Using 
previous values we thus report a luminosity of log(L/LQ) = 
0.984 + 0.04. 

Ra pid rotation modifies the location of stars in the HR di¬ 
agram. iMichel et all (I 1999 I) proposed a method to determine 
the effect of rotation on photometric parameters. In the frame¬ 
work of S Scuti stars, this method was then further devel- 
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Table 2. Table containing the main characteristics of M = I.TOMq computed models. From left to right, Q represents the 
rotational velocity in km s“'; Teff the effective temperature (on a logarithmic scale); L/Lq the luminosity in solar luminosities (on 
a logarithmic scale); R/R© the stellar radius in solar radii; g the surface gravity in cgs (on a logarithmic scale); the age in Myr; 
v(n()) the frequency of the fundamental radial mode (in pHz) and hnally, 111 /Do represents the ratio between periods of the hrst 
overtone and the fundamental radial mode. Variables indexed with a d represent the same quantities obtained when including 
near-degeneracy effects. 



n 

Teff 

L/Lq 

R/Ro 

8 

Age 

v(no) 

no/n, 

v(n;() 

(no/n,)' 

mi 

50 

3.853 

1.045 

2.184 

3.989 

1200.0 

142.553 

0.772 

142.504 

0.772 

m2 

100 

3.851 

1.042 

2.202 

3.982 

1192.5 

140.602 

0.776 

139.988 

0.773 

m3 

150 

3.849 

1.033 

2.192 

3.986 

1144.2 

141.355 

0.782 

139.176 

0.760 

m4 

200 

3.847 

1.017 

2.181 

3.990 

1039.0 

141.711 

0.792 

141.711 

0.757 

ms 

250 

3.854 

0.991 

2.046 

4.045 

830.0 

156.566 

0.798 

147.819 

0.722 

me 

50 

3.859 

1.038 

2.112 

4.018 

1132.0 

150.315 

0.773 

150.166 

0.772 

m7 

100 

3.858 

1.032 

2.101 

4.023 

1094.9 

151.447 

0.776 

150.303 

0.770 

mg 

150 

3.858 

1.020 

2.079 

4.032 

1027.5 

153.719 

0.782 

150.688 

0.767 

mg 

200 

3.850 

1.007 

2.120 

4.015 

959.0 

148.333 

0.794 

141.738 

0.728 

mio 

250 

3.859 

0.981 

1.976 

4.076 

741.0 

165.481 

0.797 

156.385 

0.753 

mu 

50 

3.867 

1.026 

2.007 

4.062 

1021.0 

162.571 

0.773 

162.616 

0.773 

mi2 

100 

3.865 

1.021 

2.013 

4.060 

999.0 

161.760 

0.776 

162.960 

0.785 

mi3 

150 

3.866 

1.006 

1.972 

4.078 

898.3 

166.858 

0.782 

170.549 

0.799 

mi4 

200 

3.859 

0.993 

2.001 

4.065 

829.0 

162.627 

0.792 

155.581 

0.758 

mis 

250 

3.864 

0.970 

1.906 

4.107 

634.0 

175.117 

0.797 

165.244 

0.752 

mie 

50 

3.874 

1.013 

1.910 

4.106 

899.4 

175.432 

0.773 

175.445 

0.773 

mi7 

100 

3.875 

1.002 

1.879 

4.119 

823.0 

179.763 

0.776 

179.988 

0.784 

mi8 

150 

3.873 

0.990 

1.872 

4.123 

755.7 

180.692 

0.782 

182.011 

0.805 

mi9 

200 

3.867 

0.980 

1.903 

4.108 

699.0 

175.909 

0.791 

182.864 

0.822 

mzo 

250 

3.867 

0.955 

1.848 

4.134 

509.0 

183.512 

0.798 

190.824 

0.885 


Table 3. Idem Table|2for M 

= 1.75M 

0 models. 








Q. 

Teff 

L/L© 

R/Ro 

8 

Age 

v(no) 

no/n, 

v(n^) 

(no/n,)‘ 

ni2i 

50 

3.850 

1.109 

2.384 

3.925 

1211.3 

126.616 

0.773 

126.598 

0.772 

m22 

100 

3.849 

1.105 

2.383 

3.926 

1192.5 

126.544 

0.776 

126.205 

0.771 

m23 

150 

3.849 

1.095 

2.356 

3.935 

1144.2 

128.588 

0.782 

126.980 

0.757 

m24 

200 

3.849 

1.077 

2.311 

3.952 

1040.0 

132.064 

0.792 

134.845 

0.765 

ni25 

250 

3.848 

1.046 

2.239 

3.980 

820.4 

137.080 

0.807 

147.028 

0.954 

m26 

50 

3.859 

1.101 

2.271 

3.967 

1133.8 

136.626 

0.773 

136.584 

0.773 

m27 

100 

3.859 

1.094 

2.248 

3.976 

1094.9 

138.732 

0.777 

138.042 

0.773 

m28 

150 

3.860 

1.081 

2.209 

3.992 

1027.5 

142.391 

0.783 

139.861 

0.769 

m29 

200 

3.859 

1.062 

2.167 

4.008 

911.0 

146.203 

0.792 

140.219 

0.735 

m3o 

250 

3.859 

1.037 

2.106 

4.033 

727.6 

151.842 

0.803 

151.842 

0.756 

m3i 

50 

3.870 

1.087 

2.124 

4.026 

1013.7 

151.508 

0.773 

151.667 

0.774 

m32 

100 

3.868 

1.082 

2.135 

4.021 

999.0 

150.273 

0.777 

148.854 

0.769 

m33 

150 

3.870 

1.066 

2.073 

4.047 

898.2 

157.023 

0.783 

153.231 

0.764 

m34 

200 

3.869 

1.045 

2.035 

4.063 

768.6 

161.177 

0.792 

154.113 

0.758 

mss 

250 

3.869 

1.025 

1.982 

4.086 

616.2 

167.403 

0.799 

157.393 

0.751 

mse 

50 

3.880 

1.071 

1.991 

4.082 

880.3 

167.222 

0.773 

167.240 

0.773 

m37 

100 

3.880 

1.061 

1.967 

4.093 

823.0 

170.376 

0.776 

170.733 

0.783 

m38 

150 

3.879 

1.048 

1.951 

4.100 

755.7 

172.435 

0.782 

174.374 

0.807 

m39 

200 

3.880 

1.020 

1.878 

4.132 

555.0 

182.397 

0.792 

185.953 

0.848 

nxio 

250 

3.878 

1.012 

1.878 

4.133 

497.4 

182.263 

0.796 

187.153 

0.866 
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oped bv IPerez Hernandez et alJ lll999l) . showing errors around 
200 - 300 K and 0.1 - 0.2 mag in the effective tempera- 
ture and absolute magnitude determination respectively (see 
ISuarez et al.ll2002l for recent results on 6 Scuti stars in open 
clusters). The error box shown in Fig. ^ has taken such sys¬ 
tematic errors into account and it will be the reference for our 
modelling. 

From photometric measurements, an estimate of the mass 
(1.75 + O.lMn) and the radius (1.58 Rp) of the star is given by 
IZakhozh^ ( ll979l) : IZakhozhai & ShaDarenkol(n~996l) . However, 
Altair is found to rotate quite ra pidly. Taking advantage of the 
fact that Altair is a nearby star, Richichi & Percherori (I 2 OO 2 I) 
measured its radius providing a diameter of 3 .12 mas, which 
corresp onds to a radius of 1.72 Rq. Moreover. Ivan Belle et alJ 
il20nili showed Altair as oblate. Their interferometric observa¬ 
tions report an equatorial diameter of 3.46 mas (correspond¬ 
ing to a radius of 1.9 R©) and a polar diameter of 3.037 
mas, for an axial ratio alb — 1.14 + 0.029. The same au¬ 
thors derived a projecte d rotational velocit y of vsin/ = 210 + 
13kms“'. Furthermore. iRover et al.lll2002h reported vsin/ val¬ 
ues which range which range from 190kms * to 250kms '. 
In addition to this, recent spectroscopically determined con- 
straints on Altair’s incli nation angle i have been established by 
Reiners & Rovei](l2004l) . They provide a range of / between 45° 


and 68°, yielding therefore a range of possible equatorial veloc¬ 
ities of Altair between 305 and 245kms“'. As shown in next 
sections, such velocities, representing 70-90% of the break up 
velocity, place significant limits on our ability to model the star. 



Fig. 2. Variation of the first overtone to fundamental radial pe¬ 
riod ratio as a function of the rotational velocity. For each mass 
considered, four curves are displayed, representing the results 
obtained modifying the effective temperature of models from 
Teff = 3.85 (bottom) to Tefr = 3.88 (top) in steps of 0.01 dex. 


3. Modelling the star 

3.1. Equilibrium models 


Stellar m odels have bee n computed with the evolutionary code 
CESAM (E ore 119971) . Around 2000 mesh points for model 
mesh grid (B-splines basis) as well as numerical precision is 
optimized to compute oscillations. 

The equation _ of state CEFF 

I Christensen-Dalsgaard & Daeppe3ll992h is u sed, in which, 
the C oulombian correction to the classical EFF dEggleton et al.l 
Il973h has been included. The pp chain as well as the CNO 
cycle nuclear reactions are considered, in which standard 
species from 'H to are included. The species D, ’Li and 
’Be have been set at equilibrium. Eor evolutionary stages 
considered in this work, a w eak elect r onic s creening in these 
reactions can be assumed dciavtonl Il968l). Opacity table s 
are taken from the OPAL package dlglesias & Rogersl [l9^. 
complement ed at low temperatures (T < 10^ K) by the tables 
provided by dAlexander & Eergusonll99'3) . For the atmosphere 
reconstruction, the Eddington T(t) law (grey approximation) 
is considered. A solar metallicity Z = 0.02 is used. Convective 
transport is described by the classical Mixing Length theory, 
with efficiency and core overshooting parameters set to 
a ML - Im/Hp = 1-8 and dov - Lv/ii n - 0.2, respective! 


The latter parameter is prescribed by Schaller et 

for intermediate mass stars. Hp corresponds to the local 

pressure scale-height, while 1^ and dgy represent respectively 


the mixing length and the inertial penetration distance of 
convective elements. 

Rotation effects on equilibrium models {pseudo rotating 
models) has been considere d by modifying the equations 
dKippenhahn & Weigerj[l990t) to include the spherically sym¬ 
metric contribution of the centrifugal acceleration by means of 
an effective gravity geff - g-^c{r), where g corresponds to the 
local gravity, and .!7lc(r) represents the radial component of the 
centrifugal acceleration. During evolution, models are assumed 
to rotate as a rigid body, and their total angular momentum is 
conserved. 

In order to cover the range of vsin/ given in Sect. |2] a set 
of rotational velocities of 50, 100, 150, 200 and 250 km s ' has 
been considered. The location in the HR diagram of models 
considered in this work is given in Eig. [0 A wide range of 
masses and rotational velocities is delimited by the photometric 
error box. To illustrate this, a few representative evolutionary 
tracks are also displayed. Characteristics of computed equilib¬ 
rium models (filled circles) are given in Tables |2l|3land|4] for 
models of 1.70 M©, 1.75 M© and 1.76 M© respectively. 

3.2. The oscillation computations 

Theoretical oscillation spectra are computed from the equi¬ 
librium models described in the p revious section. Eor this 
purpose the oscillation code Filou dTran Minh & T^ll99.5t 
ISuarezl2()o3) is used. This code, based on a perturbative anal- 
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Table 4. Idem Table|2for M - I .76Mq models. 



£1 

Teff 

L/L© 

R/Ro 

8 

Age 

v(no) 

Ho/n, 

yOfd) 

(no/n,)" 

m4i 

50 

3.851 

1.121 

2.410 

3.918 

1200 

124.959 

0.773 

124.943 

0.772 

0142 

100 

3.848 

1.117 

2.429 

3.911 

1192.5 

123.249 

0.776 

122.984 

0.771 

m43 

150 

3.849 

1.107 

2.394 

3.924 

1144.2 

125.873 

0.782 

124.407 

0.756 

ni44 

200 

3.847 

1.092 

2.375 

3.931 

1039.0 

126.768 

0.793 

131.618 

0.874 

m45 

250 

3.854 

1.062 

2.224 

3.988 

830.0 

139.769 

0.803 

144.078 

0.903 

nUe 

50 

3.859 

1.113 

2.306 

3.957 

1132.0 

133.931 

0.773 

133.896 

0.773 

0147 

100 

3.859 

1.106 

2.285 

3.965 

1094.9 

135.737 

0.777 

135.124 

0.773 

11148 

150 

3.860 

1.094 

2.238 

3.983 

1027.5 

139.989 

0.783 

137.570 

0.769 

11149 

200 

3.857 

1.083 

2.248 

3.979 

959.0 

138.644 

0.791 

133.851 

0.736 

11150 

250 

3.861 

1.051 

2.123 

4.029 

741.0 

150.734 

0.801 

150.734 

0.758 

msi 

50 

3.870 

1.100 

2.161 

4.013 

1021.0 

148.028 

0.773 

147.652 

0.771 

m52 

100 

3.868 

1.095 

2.165 

4.012 

999.0 

147.557 

0.777 

146.313 

0.770 

11153 

150 

3.870 

1.078 

2.096 

4.040 

898.3 

154.935 

0.783 

151.293 

0.764 

0154 

200 

3.865 

1.066 

2.123 

4.029 

829.0 

151.476 

0.792 

145.163 

0.759 

0155 

250 

3.868 

1.039 

2.024 

4.070 

634.0 

162.494 

0.801 

152.321 

0.750 

m56 

50 

3.879 

1.085 

2.031 

4.067 

899.4 

162.756 

0.773 

162.782 

0.773 

OI 57 

100 

3.881 

1.073 

1.988 

4.086 

823.0 

168.121 

0.777 

168.527 

0.783 

11158 

150 

3.880 

1.060 

1.968 

4.095 

755.7 

170.675 

0.783 

172.777 

0.808 

OI 59 

200 

3.875 

1.051 

1.995 

4.083 

699.0 

166.886 

0.791 

159.662 

0.756 

moo 

250 

3.873 

1.021 

1.944 

4.105 

509.0 

172.969 

0.802 

161.557 

0.749 


ysis, provides adiabatic oscillations corrected for the effects of 
rotation up to second order (centrifugal and Coriolis forces). 

Furthermore, for moderate-high rotational velocities, the 
effects of near d egeneracy are expected to be signihcant 
llSouh et alJl998ll . Two or more modes, close in frequency, are 
rendered degenerate by rotation under certain conditions, cor¬ 
responding to selection rules. In particular these rules select 
modes w ith the same azim uthal order m and degrees £ differ¬ 
ing by 2 dSoiifi etalJll998l> . If we consider two generic modes 
ai = (n,£,m) and q '2 = {n',£',m') under the aforementioned 
conditions, near degeneracy occurs for Icr^j -cr^J < ctq, where 
CTa, and 0 - 0,2 represent the eigenfrequency associated to modes 
a I and a '2 resp ectively, and crp re presents the stellar rotational 
frequency tsee lnoiinil et al .1200(1 for more details). 

4. Comparison between theory and observations 

High-amplitude 6 Scuti stars (HADS) display V amplitudes in 
excess of 0.3 mag and generally oscillate in radial modes. In 
contrast, lower-amplitude members of the class present com¬ 
plex spectra, typically showing non-radial modes. 

The amplitude of the observed main frequency is below 
0.5 ppt which is several times smaller than typically detected 
in non-HADS 6 Scuti stars from ground-based observations. 
As shown in paper I, using the classical period-luminosity re¬ 
lation with the fundamental parameters giv en in Sect. E| and 
assuming the value of Q - 0.033 d given bv lBrege3lll979li for 
6 Scuti stars, the frequency of the fundamental radial mode is 
predicted to be 15.433 d“', suggestively close to the observed 
/i = 15.768 d '. For the remaining modes (Tabled there is no 



-i.te -1.14 -1.12 

log rio (d) 



-1.12 -1.10 -1.08 -1.06 -1.04 

log Ho (d) 



-1.12 -1.10 -1.08 -1.06 
log Ylo (d) 


Fig. 3. Hi/Ho versus logHo. Three panels are displayed corre¬ 
sponding to the three masses considered in this work. For each 
panel, different lines represent the results obtained varying the 
rotational velocity from Q = 50 km s ' (bottom, solid line) to 
Q = 250kms“' (top, dash-dotted line). 
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n (km/s) 

Fig. 4. Same as Fig. |2lbut considering near degeneracy. Te^i and Te ,2 labels represent effective temperatures expressed in a 
logarithmic scale in K. The vertical thick solid line indicates the limit for which Fti /Do varies linearly with Q. in presence of 
near degeneracy (see details in text). The horizontal dashed thick line represents the observed fli/Ilo. Finally, the shaded area 
represents the range of rotational velocities Q. (assuming i = 90 deg). 


observational evidence to identify them as radial. Nevertheless, 
the second dominant mode, / 2 , will be considered here as the 
first overtone. 


The present study is divided into two parts. In Sect. 14.11 
we consider only the two dominant modes (ie. f\ and / 2 , cf. 
Table The observations are then compared with models 
through period-ratio vs. pe riod diagrams, from now on called 


Peter s en dia grams (see e.g. [Petersen & Christensen-Dalsgaard 
1 19961 Il999l) . Then, in Sect. lO we use the 5 most dominant 
observed modes to find the best fit to the theoretical models. 


4.1. Considering two radiai modes 


The well known Petersen diagrams show the variation of 
Hi/Do ratios with log Ho, where Ho represents the period of 
the fundamental radial mode, and IIi the period of the first 
overtone. They are quite useful for constraining the mass and 
metallicity of models given the observed ratio between the fun¬ 
damental radial mode and the first overtone. However these di- 
agra ms do not consider the effect of rotation on o scillations 
(see Petersen (^Christensen-Dalsgaard|[l996l 1 19991 for more 
details). Problems arise when rapid rotation is taken into ac¬ 


count; it introduces variations in the period ratios even for ra¬ 
dial modes. This, combined with the sensitivity of Petersen di¬ 
agrams to the metallicity, renders the present work somewhat 
limited. A detailed investig ation of such combined effects is 
currently a work in progress ISuarez & Garridol2005l) . In order 
to cover the range of effective temperature (from Teff = 3.85 
to Teff = 3.88), luminosity and rotational velocity, three sets 
of 20 models are computed, where rotational velocity varies 
from Q = 50 to Q = 250kms '. Characteristics of the three 
sets are given in Tables |2] 0 and 0 corresponding to sets of 
I.TOMq, 1.75Mo and 1.76Mo models respectively. Adiabatic 
oscillations are then computed considering radial and non- 
radial modes. In order to analyse the Petersen diagrams cor¬ 
responding to these models we calculate the Hf/Ho ratios. In 
Fig. El the Of/Ho ratios are displayed as a function of the ro¬ 
tational velocity Q. For each mass (line style) there are four 
curves corresponding to four different effective temperatures 
(from bottom to top). It can be seen that all ratios increase pro¬ 
portionally with the rotational velocity. From Q = 50 km s ' 
to Q = 200 km s ' all models behave identically following a 
^MdSi-linear relation. For low rotational velocities, Hi/Ho ra¬ 
tios approach the canonical 0.77 obtained in absence of rota- 
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3.aao 3.a75 3.a70 s.aes s.aeo 3.a.55 3.a50 3.aao 3.a75 3.a70 s.aes s.aeo 3.a55 3.a50 

logT^„ (K) logT,„ (K) 

Fig. 5. Variation of fundamental mode frequencies v(no) (in yuHz) with effective temperature for models with a rotational velocity 
of Q = 200 km s’, representative of Altair. Left panel corresponds to frequencies not corrected for the effect of near degeneracy, 
and right panel to those corrected. Finally, the horizontal dashed line represents the frequency of the observed fundamental mode 

rob.s 

JQ • 


tion. As expected, dispersion of curves is not signihcant, since 
the masses of the models are quite similar, for a given rotational 
velocity. However, for Q. > 200kms“’, dispersion of curves 
also increases as a consequence of errors in the perturbative 
method used for computing oscillations. With rotational veloc¬ 
ities higher than ~ 200 km s ’, results should be interpreted 
carefully, since the se cond order pert urbation theory consid¬ 
ered here may fail (see Lignieresl200l[ for more details). With 
a projected velocity of vsini = 190-230kms^’, Altair is at the 
limit of validity. 


The impact of these results on Petersen diagrams are shown 
in Fig. 13 In this diagram the variation of Hi/Ho is displayed 
versus log Hq for the three masses studied. As in Fig. |2] note 
the trend toward the standard value of Hi /Ho = 0.77 as the ro¬ 
tational velocity decreases. Considering the Hi /Ho dependence 
in isolation, a canonical value Hi /no(Q) can be established for 
a given mass. For the range of masses studied here, the analysis 
of three panels of Fig.|3shows a very low dependence on mass. 


Up to this point, the effect of near degeneracy has not been 
taken into account. Near degeneracy occurs systematically for 
close modes (in frequency) following certain selection rules 
(see Sect. m. It increases the asymmetry of multiplets and 
thereby the behaviour of modes. The higher the value of the ro¬ 
tational velocity, the higher the importance of near degeneracy. 
In the present case, for the range of rotational velocities con¬ 
sidered, near degeneracy cannot be neglected. As can be seen 
in Fig. |4] not all models present the same behaviour with Q. 
For rotational velocities up to 150kms ’, a double behaviour 


is shown, one for lower effective temperatures (FIj /Hq < 0.77), 
and one for higher effective temperatures (H^/Ily > 0.77), 
both remaining linear. This shows the dependence of Hi /no(Q) 
with the evolutionary stage of the star. However, for higher 
rotational velocities, particularly on the right side of the ver¬ 
tical limit line, everything becomes confusing. At these ro¬ 
tational velocities, the global effect of rotation up to second 
order (which includes asymmetries and near degeneracy) on 
multiplets complicates the use of H^/nQ(Q) predictions on a 
Petersen diagram. In particular, radial modes are affected by 
rotation through the distortion of the star (and thereby its prop¬ 
agation cavity), and thro ugh near degenera cy effects, coupling 
them with { - 2 modes (ISoiifi et alJll998l) . In addition, other 
factors such as the evolutionary stage and the metallicity must 
also be taken into account, which makes the global dependence 
of Petersen diagrams on rotation rather complex. 


In Fig.|3(left and right panels), the combined effect of rota¬ 
tion and evolution on the fundamental radial mode is shown for 
representative models with a rotational velocity of 200kms“’. 
In particular, for the 1.76 Mq model, the stellar radius is ap¬ 
proximately 10 ’ Rq larger than those of 1.70 and 1.75 Mq, 
which difference is of the order of 10“^ R©. As a result, a clear 
difference of behaviour between those models is observed. 
Such a difference depends not only on the mass, but also on 
the evolutionary s tage, the inital rotati onal velocity (at ZAMS), 
metallicity, etc. In lSiiarez et al.lil200.5l) . for a certain small range 
of masses, a similar non-linear behaviour is found for predicted 
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Table 5. Table containing a proposed mode identification for the observed frequencies given in Table [2 using the oscillations 
computed from different models. Values in parentheses correspond to the mode identihcation numbers (n, i, m). The hrst hve 
models corresponds to those chosen by their proximity to the fundamental radial mode. The last hve, corresponds to those 
chosen by minimum cr^ value, v, represent the theoretical frequencies closest to the observed ones. 




< 

V2 

2 

V3 

^3 

V4 


Vs 


(T^ 


mi8 

(-1,2,0) 

(1,0,0) 

(1,2,2) 

(1,2,2) 

(2,2,-l) 

(2,2,-l) 

(-1,2,0) 

(-1,2,0) 

(-1,2,-!) 

(-1,2,-!) 

16.937 

5.698 

mi9 

(-2,1,0) 

(1,0,0) 

(3,2,2) 

(3,2,2) 

(4,0,0) 

(4,0,0) 

(-3,2,-2) 

(-1,1,0) 

(-1,1,-!) 

(-3,2,-2) 

12.531 

4.056 

m2o 

(1,0,0) 

(-1,2,0) 

(1,2,2) 

(1,2,2) 

(2,2,-l) 

(2,2,-l) 

(-1,2,0) 

(-1,2,0) 

(-1,2,-!) 

(-1,2,-!) 

4.575 

4.544 

m39 

(1,0,0) 

(1,1,0) 

(2,2,1) 

(2,2,1) 

(6,2,2) 

(2,2,0) 

(-2,2,-l) 

(-2,2,-l) 

(-2,2,-l) 

(-2,2,-1) 

21.234 

7.600 

m4o 

(1,0,0) 

(-l,2,-2) 

(1,2,0) 

(1,2,0) 

(3,1,1) 

(3,1,1) 

(-l,2,-2) 

(-l,2,-2) 

(1,1,1) 

(1,1,1) 

0.763 

1.242 

mio 

(-1,1,0) 

(-2,2,0) 

(3,2,2) 

(3,2,2) 

(4,0,0) 

(4,0,0) 

(-3,2,-2) 

(-1,1,0) 

(-1,1,-!) 

(-3,2,-2) 

1.194 

0.820 

mi3 

(-1,2,0) 

(-1,2,0) 

(1,2,2) 

(1,2,2) 

(2,2,-l) 

(2,2,-l) 

(-1,2,0) 

(-1,2,0) 

(-1,2,-!) 

(-1,2,-!) 

0.903 

0.903 

m23 

(-l,2,-2) 

(-l,2,-2) 

(3,2,2) 

(3,2,2) 

(4,1,-1) 

(4,1,-1) 

(0,2,0) 

(0,2,0) 

(0,2,0) 

(0,2,0) 

0.835 

0.835 

m34 

(-1,2,0) 

(1,1,0) 

(2,2,1) 

(2,2,1) 

(6,2,2) 

(2,2,0) 

(-2,2,-1) 

(-2,2,-l) 

(-2,2,-l) 

(-2,2,-1) 

1.698 

1.228 

m58 

(-l,2,-2) 

(-l,2,-2) 

(2,1,1) 

(2,1,1) 

(3,1,-1) 

(3,1,-1) 

(-1,1,1) 

(-1,1,1) 

(1,1,1) 

(1,1,1) 

0.639 

0.940 


unstable mode ranges. The results of this work may provide a 
clue to understand this unsolved question. 

The reader should notice that at these velocities, when fre¬ 
quencies are corrected for near degeneracy effects, their vari¬ 
ations with Teff (and thereby with evolution) are more rapid 
(right panel) than was the case for uncorrected ones. Moreover, 
a shift to higher frequencies is observed when correcting for 
near degeneracy, which favors the selection of lower mass ob¬ 
jects. For low temperatures (i.e. for more evolved models), sta¬ 
bilization is found for low v(nQ) frequencies, which can be ex¬ 
plained by the joint action of both near degeneracy and evolu¬ 
tion effects. 

In the present work, following the prescription of 
iGoupil et alJ (l2000l) . modes are near degenerate when their 
proximity in frequency is less or equal to the rotational fre¬ 
quency of the stellar model (Iv,- - Vj\ < ctq). Considering all 
possibilities, 5 models identify the observed /i as the funda¬ 
mental radial mode; mig, mig, m 2 o (1.70Mq), mgg and m^o 
(1.75 Mq). For last three models (20, 39 and 40), this happens 
without considering near degeneracy effects. For mig and mig, 
the theoretical fundamental mode approaches /i when consid¬ 
ering near degeneracy. 

In Table 13 radial and non radial identihcations (discussed 
in the next section) are presented for selected models. The hrst 
set of hve models corresponds to those selected by their identi¬ 
hcation of /i as the radial fundamental mode. As can be seen, 
no identihcation is possible when trying to ht the whole set 
of observed frequencies. However /g is identihed as the third 
overtone by the rapid rotating model (Q = 200 km s ') mig. 
Uncertainties in the observed mass and metallicity are also an 
important source of error in determining the correct equilib¬ 
rium model. Thus, the fact of obtaining fundamental modes 
with frequencies lower than 180/iHz for most of the mod¬ 
els could be explained by an erroneous position of the pho¬ 
tometric box on the HR diagram. In fact, the lower the mass 
of model used (always within the errors), the higher the v(no) 
value. However, the location in the HR diagram of models with 


masses lower than 1.70 Mq (with the same metallicity) is not 
representative of the Altair observations. 

At this stage, there are two crucial aspects to hx. On 
one hand, it is necessary to determine the physical conditions 
which enforce degeneracy between mode frequencies. A phys¬ 
ically selective near-degeneracy could explain the behaviour of 
Hi/Ho for very high rotating stars. On the other hand, for high 
rotational velocities, third order effects of rotation are presum¬ 
ably important. 

4.2. Non-radial modes. 

As neither observational nor theoretical evidence supporting 
a radial identihcation of the observed modes fj exists, in this 
part of the work, we carry out an analysis generalised to non- 
radial modes. As we did for radial oscillations (see Sect, 
here we compute non-radial adiabatic oscillations for the mod¬ 
els of Tables |2j |3 and IS Assuming the 7 observed modes of 
Altair pulsate with ^ < 3, a possible mode identihcation is pro¬ 
posed in Table |3 In order to avoid confusion, only the hrst 5 
observed modes (Tabled) with larger amplitudes, will be con¬ 
sidered. This constitutes a rough identihcation based only on 
the proximity between observed and theoretical mode frequen¬ 
cies for each model. That is, no additional information about 
the £ and m values is used. 

In order to obtain an estimate of the quality of hts (iden¬ 
tihcation of the whole set of observed frequencies), the mean 
square error function is used. The lower its value, the better the 
ht for the free parameters considered. For each model 

is calculated, where foj and v, ,■ represent the observed and the¬ 
oretical frequencies respectively. The total number of observed 
frequencies is represented by A = 5. Calculations are made 
hxing the metallicity Z, the overshooting parameter dov and the 
mixing length parameter auh (see Sect. 13. H . On the other side. 
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Fig. 6 . HR diagram showing selected models given in Tabled Labels correspond to model numbers. Square symbols represent 
models selected by the proximity to the fundamental radial mode; filled squares are used for those taking into account near 
degeneracy effects. Filled circles represent models selected by the minimum cr^ value. As a reference, observations (the cross 
symbol) and two evolutionary tracks are depicted; One corresponding to a 1.70Mo and Qi = 150kms“' (solid line), and the 
other one, corresponding to a 1.75 Mq and O 2 = 200 km s ' (dashed line). For more details see the text. 


the mass M, the rotational velocity Q. and the evolutionary stage 
Xc have been considered as free parameters. In Tabled the last 
two columns list cr^(M, Q, Xc) calculated for each model, both 
without near degeneracy (cr^) and including it (cr^). When con¬ 
sidering the whole set of models computed, the inclusion of 
near degeneracy increases the quality of fits by roughly 30%. 
Table 13 is divided in two parts: the first 5 models have been 
chosen as identifying fi as the fundamental radial mode. The 
second set of models correspond to those with minimizing cr^. 
Analysing the results obtained for the 14 selected models it can 
be seen that no ^ = 3 identifications are found and most of 
identifications correspond to { — 2 ( 60%) and £ - 1 ( 32%) 
modes. As happened for the first set of models, /a is identified 
as the third radial overtone by models miQ. 

On the other hand, since no specific clues for the degree £ 
and the azimuthal order m are given, other characteristics of the 
observed spectrum must be investigated. In particular, it is quite 
reasonable to consider some of the observed frequencies as be¬ 
longing to one or more rotational multiplets. Although no com¬ 
plete multiplets are found, several sets of two multiplet mem¬ 
bers are obtained. Specifically, the observed frequencies and 
/s are identified as members of ^ = 2 multiplets by models 
mi 3 , mis, ni 2 o, and mio, and as members of £ - \ triplets 
by models mio and mig. 

Figure |3 shows the location on a HR diagram of selected 
models given in Table |3 Notice that around 80% of them lie 


in effective temperature and luminosity ranges of log Teif = 
[3.88,3.85] and logL/Lo = [0.96,1.09] respectively. In the 
area delimited by these models (shaded surface), models with 
similar characteristics are found. As can be seen, the shaded 
area is located in the central part of the error box. Except mod¬ 
els m 39 , m 4 o and mss, all models are situated toward colder 
and more luminous locations in the box. This is in agreement 
with the expected effect of rotation on fundamental parameters 
(see Sect. 13 . On the other hand, positions of models selected 
by the proximity to the fundamental radial mode (squares) can 
be connected by a straight line. This is an iso-p line', with 
p ^ 5.23 gcm^, and which basically explains the similar fre¬ 
quencies of their fundamental radial mode. Considering near¬ 
degeneracy, the identification of the fundamental radial mode 
and the lowest cr^ value, our best model is mig, which is at the 
lower limit of the range of Altair’s observed vsin/. 

Nevertheless, in order to further constrain our models rep¬ 
resentative of Altair, it would be necessary to obtain additional 
information on the mode degree and/or azimuthal order of ob¬ 
servations. In this context, spectroscopic analysis may provide 
information about £, m and the angle of inclination of the star, 
as nonradi al pulsations generat e Doppler shifts and line profile 
variations jAerts & FvenCOOfli . Furthermore, multicolor pho- 


' Line which connects models with the same stellar mean density. 
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tomet ry may also provide information about £ dGarrido et alJ 
1990ll . 

5. Conclusions 

In the present paper a theoretical analysis of frequencies of 
HR6534 (Altair) was presented, where rapid rotation has been 
properly taken into account in the modelling. The analysis was 
separated in two parts; 1) considering the observed modes /i 
and /2 corresponding to the fundamental radial mode and the 
first overtone and models were analysed through the Petersen 
diagrams, and 2) a preliminary modal identification was pro¬ 
posed by considering radial and non-radial oscillations. 

Firstly, in the context of radial modes, we studied the iso¬ 
lated effect of rotation on Petersen diagrams. For the different 
rotational velocities considered, the shape of 111 /!!()(□) leads 
to a limit of validity of the perturbation theory (up to second 
order) used at around Q = 200kms '. This limit is mainly 
given by the behaviour of such period ratios when near degen¬ 
eracy is considered, which visibly complicates the interpreta¬ 
tion of Petersen diagrams. Nevertheless, in this procedure, for 
rotational velocities up to 150kms“' it is found that Ili/Ilo 
is lower than 0.77 for lower effective temperatures, and recip¬ 
rocally higher than 0.77 for higher effective temperatures (in¬ 
side the photometric error box). The analysis of radial modes 
also reveals that only a few models identify the observed /i 
as the fundamental radial mode. This reduces the sample of 
models to those with masses around 1.70Mo or lower, on the 
main sequence and implies models with different metallicity 
in order to agree with photometric error box. These partial re¬ 
sults (for a given region in the HR diagram, for a given range 
of masses, evolutionary stages and a given metallicity) consti¬ 
tute a promising tool for seismic investigation, not only of 6 
Scuti stars, but in general, of multi-periodic rotating stars. A 
detailed investigation on the effect of both rotation and metal- 
licitv on Petersen diagram s will be proposed in a coming paper 
dSuarez & Garridol2005l) . 

Secondly, in the context of radial and non radial modes and 
assuming observed modes pulsating with £ = 0,1,2,3, a set of 
14 models was selected in which five of them identify the fun¬ 
damental radial order and at least six others identify two of the 
observed frequencies (/4 and /s) as members of ^ = 2 and £ - I 
multiplets. A range of masses of [1.70,1.76] Mq principally for 
a wide range of evolutionary stages on the main sequence (ages 
from 225 to 775 Myr) was obtained. Considering informa¬ 
tion of both radial (through Petersen diagrams) and non-radial 
modes, a set of representative models with rotational velocities 
larger than 150 km s ' was obtained. 

Further constraints on the models are thus necessary. Such 
constraints can be obtained by employing additional informa¬ 
tion on the mode degree £ and/or the azimuthal order m of the 
observed modes, which may be inferred from spectroscopy and 
multicolor photometry dGarrido et al ■Il990l) . Improvements on 
adiabatic oscillation computations (including third order com¬ 
putations) will constitute a coherent and very powerful tool to 
obtain seismic data from future space missions like COROT, 
EDDINGTON and MOST. 
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